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Abstract. We discuss a model of X-ray variability of active galactic nuclei (AGN). We consider multiple spots 
which originate on the surface of an accretion disk following intense irradiation by coronal flares. The spots 
move with the disk around the central black hole and eventually decay while new spots continuously emerge. We 
construct time sequences of the spectra of the spotted disk and compute the corresponding energy-dependent 
fractional variability amplitude. We explore the dependence on the disk inclination and other model parameters. 
AGN seen at higher inclination with respect to the observer, such as Seyfert 2 galaxies, are expected to have 
fractional variability amplitude of the direct emission by a factor of a few higher than objects seen face on, such 
as the Seyfert Is. 
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1. Introduction 

Broad band spectra of active galactic nuclei show the 
presence of several components. The two principle con- 
tributions are: (i) the Big Blue Bump extending from 
optical/UV band to (sometimes) soft X-ray band, and 
(ii) the hard X-ray power law. The Big Blue Bump emis- 
sion is conventionally interpreted as originating from an 
accretion disk (Czerny & Elvis 1987; Koratkar & Blaes 
1999; however, see e.g. Collin et al. 2002 for discussion 
of problems encountered by this scheme). On the other 
hand, the nature and geometry of the region responsi- 
ble for the hard X-ray emission is still under discussion. 
Several models were proposed, the most popular ideas 
being a hot extended corona overlaying a relatively cold 
disk (e.g. Czerny & Elvis 1987, Haardt & Maraschi 1991, 
Rozanska & Czerny 2000, Liu et al. 2002, Merloni 2003), 
an inner hot flow (e.g. Ichimaru 1977, Narayan & Yi 1994, 
Narayan et al. 2002) the lamp-post model (where a point- 
like X-ray source is located at a specified height above the 
disk, e.g. Henri & Pelletier 1991, Malzac et al. 1998), and 
the model of multiple hot flares produced via magnetic 
field reconnections (Galeev et al. 1979 and later papers). 
The field was reviewed, for example, by Leighly (1999), 
Collin (2001), and Poutanen (1998) and Done et al. (2002) 
(in the context of X-ray binaries). Notice that these sce- 
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narios are not completely disparate and there may be a 
certain overlap between them. In all models the primary 
source of X-rays should be strongly variable. 

Indeed, AGN are variable in the X-ray band (e.g. 
Lawrence et al. 1987; Taylor et al. 1993) but the origin 
of variability cannot be examined directly, because the 
relevant central regions still remain unresolved to direct 
imaging. However, studies of the X-ray spectral variabil- 
ity offer a direct insight into the structure of accretion 
flows onto a central black hole, which is assumed to power 
AGN and determine their spectra. 

In the flare model, sudden dissipation occurs in very 
localized regions above the disk surface (coronal loops, 
in analogy with the solar corona). In this case the local 
irradiation flux can be orders of magnitude higher then 
the steady level of the disk flux itself, but the irradiation 
does not last very long and only a small fraction of the 
disk surface is irradiated at every moment. 

Such arrangement of X-ray emitting region has been 
suggested for the flrst time by Galeev et al. (1979), 
and it was subsequently developed in many papers (e.g. 
Abramowicz et al. 1991; Haardt et al. 1994; van Oss et 
al. 1993; Poutanen & Fabian 1999; Zycki 2002). The de- 
velopment of a flare leads to a burst of 'primary emission' 
as well as to the formation of a hot spot underlying the 
flare where roughly half of the X-ray flux is reprocessed 
by the disk. Irradiation of the disk surface in hydrostatic 
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equilibrium leads to the formation of strongly stratified 
medium - a hot fully ionized skin covering a cooler, more 
neutral zone - and so the X-ray spectrum resulting from 
reprocessing should contain signatures that are charac- 
teristic for multi-temperature gas (Nayakshin et al. 2000; 
Ballantyne et al. 2001; Rozahska et al. 2002). The irradi- 
ating flux is locally very large, exceeding considerably the 
stationary energy flux which is dissipated inside the disk 
(e.g. Nayakshin 2000; Ballantyne et al. 2001). The result- 
ing ionized skin is relatively optically thick (Collin et al. 
2003). 

The flare/spot model is a possible (although not 
unique) explanation of narrow emission features, which 
have been reported in ~ 5-6 keV X-ray spectra of sev- 
eral AGNs and interpreted in terms of in terms of local- 
ized iron-line emission (Turner et al. 2002, 2004; Guainazzi 
2003; Yaqoob et al. 2003; Dovciak et al. 2004). 

In the present paper we test the flare/spot model by 
analyzing its predictions for the fractional variability am- 
plitude in the X-ray band. We model the local spot/flare 
spectrum as a sum of flare (primary) emission of a power 
law shape, and a spot (reflected) emission. The spot emis- 
sion is determined by using the coupled titan/noar 
codes (Dumont et al. 2000) to solve the radiative transfer 
and the code of Rozahska et al. (1999) to calculate the hy- 
drostatic equilibrium. We assume random distribution of 
spots and flares across the disk and we account for their 
motion during the time-integrated observation. We con- 
sider both non-rotating (Schwarzschild) and rapidly ro- 
tating (Kerr) black holes, and we apply general relativity 
corrections using the KY code of Dovciak et al. (2003). 

In our approach, a sequence of solutions is specified 
by fixing statistical properties which the flare/spot dis- 
tribution is required to obey. Given a particular solution, 
the corresponding energy-dependent fractional variability 
amplitude Fvar is calculated. The model is described in 
Section [21 Results are given in Section 13 The relevance 
of the model for explaining spectrum and variability of a 
typical Seyfert 1 galaxy is discussed in Section^ 

2. The model 

2.1. Intrinsic spectrum of the spot 

We model the X-ray spectrum from an irradiated accretion 
disk around a massive black hole, assuming a power law 
profile of the flare primary emission and computing the 
radiation spectrum from a hot spot underlying the com- 
pact flare. The spot emission goes directly to the observer 
without any further secondary Comptonization within the 
flare (we neglect size of the flare itself) and we adopt plane 
parallel symmetry for local computations. 

Reprocessing of the incident flux is computed with the 
code TITAN (Dumont et al. 2000). The radiative transfer 
is solved both for the continuum and for lines, because the 
escape probability approach does not give reliable results 
for the lines (Dumont et al. 2003). The Comptonization 



is taken into account by coupling with the Monte Carlo 
code NOAR (also Dumont et al. 2000). 

The vertical structure of the underlying disk is cal- 
culated with the code of Rozahska et al. (1999). As our 
principal model, we consider an irradiated disk in hydro- 
static equilibrium, which requires iterations between the 
disk structure and the radiation transfer computations. 
These iterations are carried out as described by Rozahska 
et al. (2002). The computations of equilibrium are ex- 
tremely time-consuming because of coupling between the 
equations of radiation transfer for the lines and for contin- 
uum. Therefore, in the present paper we performed com- 
putations at a certain representative radius in the disk 
{R = 18i?g in terms of gravitational radii, Rg=GM/c^) 
with the aim of obtaining a typical reflection spectrum, 
F\oc{E)- Local reflection spectra at other radii were ob- 
tained with the help of scaling. 

2.2. Assumptions about the flare - spot connection 

We intend to represent a typical observation of an AGN, 
and therefore we choose as model parameters preferen- 
tially those properties which are connected as close as 
possible to observations. We define a very simple model 
(in fact, too simple to explain the origin of power spec- 
tra themselves). The model takes into account the basic 
observational properties, such as mean luminosity and lu- 
minosity dispersion. This approach is complementary to 
the papers devoted to the development of an isolated flare 
(e.g. Bottcher et al. 2003). 

Fundamental parameters of the model are the mass of 
the black hole, M, and the mean X-ray luminosity, Lx- 
This X-ray luminosity is assumed to originate from n^can 
flares per second, on average. The actual number of flares 
fluctuates at any given moment around this mean value 
according to Poisson's distribution. Flares are assumed to 
corotate with the underlying disk. 

A fraction of radiation from each flare reaches the ob- 
server as the direct component. On the other hand, a part 
of flare photons is reprocessed by the disk surface, pro- 
ducing the reflection component. 

In the model, a flare occurs at height h above the 
disc, it irradiates the disk surface and creates a hot spot. 
Irradiation is the largest and strictly perpendicular just 
below the flare (in the centre of the spot where it first 
appears). Away from the centre, the irradiation slowly de- 
creases over distance up to Rx~h from the central point 
(with inclination angle of the incident X-rays increasing) 
and it drops down fast further out. We start by simpli- 
fying our analysis and assuming that the irradiation does 
not show any gradient across the spot radius up to i?x- 
Then it drops to zero sharply. Linked with the spatial 
dependence of illumination across the spot surface there 
should be also a corresponding change in the incident an- 
gle. However, this incident angle is about 60° at most of 
the spot surface so we assume illumination to be isotropic. 
In this way, we consider the reprocessed radiation as com- 
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ing from uniformly radiating circular spots. Given the in- 
cident radiation flux of a single flare, i^inc, and its intrinsic 
spectrum Fioc (E) , the reflection is calculated through sim- 
ple scaling, 



^ inc 
^loc ' 



(1) 



where Fio^ is the radiation flux assumed in titan/noar 
computations. This approach neglects aberration and light 
bending in the local irradiation event, which is roughly 
correct if the flare height is small in comparison with the 
radius in the disk, fn this manner, we can connect the flare 
luminosity, Lj, the incident radiation flux and the size of 
the spot: 



(2) 



We further assume that all flares occur at the same 
height above the disk surface, h. This is partially moti- 
vated by the fact that the disk thickness does not de- 
pend on radius in standard, radiation pressure dominated 
Shakura-Sunyaev model, but mainly by simplicity of the 
assumption. Therefore, all spots are of the same size (i?x)- 
Finally, we assume that flares are generated above the disk 
according to uniform random distribution between inner 
radius, i?in, and outer radius, i?out- Gravitational energy 
available in the disk depends on radius, and so we cannot 
expect all flares to be identical. In the present paper we 
adopt Li as a function of flare location, in the form 



Li{Ri) (X ( 1^ 



(3) 



This formula implies that the incident radiation flux scales 
with distance of the flare from the centre as a power law. 
We denote the corresponding power-law index /Jrad and 
we normalize the flux amplitude to Fq at the inner rim. 
Therefore, 



Fo 



(4) 



Spots are instantaneously created by the flare occur- 
rence, and hence one can think of them as being attached 
to flares and orbiting with them at the local Keplerian 
speed around the central body. We assume that duration 
of all flares is the same and equal to info, and we deter- 
mine the average spectrum as seen by the observer after 
integration over the period Tobs. 

In conclusion of this section, we can summarize all nine 
input parameters of the model: M, nmean, Rim Rout, Rx, 
Fq, Aad: ^lifc: and Tobs- These parameters fully define the 
flare distribution, including the total mean luminosity of 
the source, Lx- In practice it may be more convenient to 
use Lx as a free parameter instead of, for example, Rx- 

Prediction of the observed spectrum involves two ad- 
ditional parameters which are inherent to any model in- 
volving an accretion disk around a rotating black hole: 
inclination angle of the observer, i, and dimension-less an- 
gular momentum of the black hole, a. These parameters 



span the range 0° < i < 90° {i = 0° the disk axis, while 
i — 90° is the disk plane), and 0<a<I(a = refers 
to a non-rotating black hole, while a — 1 corresponds to 
a Kerr black hole with maximum rotation). 

2.3. Properties of the flare distribution 

In previous section we introduced several assumptions 
about the flare/spot connection. In the present section we 
use these assumptions to determine analytical relations 
describing secondary properties of the flare distribution. 
Such formulae are very handy for rough estimates and for 
discussion of properties of real sources. Since we assume 
that the spot radius is small in comparison with the over- 
all size of the disk, also general relativity corrections are 
assumed to be negligible for the spot local properties.^ 

We first derive the mean luminosity, Lx- Uniform dis- 
tribution of flares over the disk surface means that prob- 
ability for a flare to occur at a given radial and azimuthal 
position, {Ri, 0i), is equal to 



p(RiAi)=p{Ri)p{c^i) 



2Ri 



??2 
^out 



Rt 



1 

2^' 



(5) 



One can associate the average total luminosity of a source 
with an exemplary distribution of flares. 



L = ^Li(i?i,0i). 

i=l 



(6) 



The mean (expected) value of this quantity is determined 
as integral over the probability distributions. 



L 



X = P{n; nmean 

n=l 



Ox 



-R.OU 



2tt 



27r 



p{Rl)...p{Rn) 



(7) 



p(0i)...p(0„) di?i... di?„ d(/)i... d(/)„ ^ (j)i), 

i=l 

where P(n; rimcan) denotes Poisson's distribution around 
"-mean- The multiple integrals can be performed easily, be- 
cause there is no dependence of the integrand on 0i. All 
integrals over Ri are identical (there are n of them) and 
can be performed independently of each other. Summation 
over the Poisson distribution gives 



n=l 



(8) 



and the expression for the mean luminosity reduces to 
a single integral over the adopted radial distribution of 
luminosity, 

(•-Rout 

ix = ?^moan / p(i?i)Li di?i. (9) 
7Ri„ 



^ It will be interesting to relax this assumption in future 
work. Also, it should be emphasized that we do not ignore 
general relativity effects acting on reprocessed photons of the 
spot (see Section . 
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This integral can be calculated analytically. Taking into 
account equations ©-lEll and denoting ^=i?in/i?out, we 
obtain 



flout 



(10) 



Using this formula, we can express the spot radius as a 
function of Lx and other parameters involved: 



Rx 



Lx (/3rad-2)(l-C') 



1/2 



(11) 



At any given moment, the mean covering factor of the disk 
surface with spots is 



^out 



Rl' 



(12) 



The disk surface is completely covered with reprocessing 
spots if Cmoan IS closc to (or greater than) unity. 

The mean number of flares during observation of du- 
ration Tobs is given by 



life 



1 



(13) 



Obviously, if the observation is much shorter than the flare 
duration then A^moan is determined just by the mean num- 
ber of flares. On the other hand, Amoan linearly increases 
with time in the opposite limit of a very long observation. 

Since some of the flares existing during the observation 
were born before the beginning of the observation or they 
will disappear after the end of observation, the mean ob- 
serving time of any single flare is clearly shorter than the 
flare duration time. We therefore introduce effective life- 
time, ij^^p, which is a combination of the two timescales: 



'-life — 



^obs ^life 
^obs ~t~ ^life 



(14) 



The model specifies only the mean number of flares, 
but in each specific observation we expect various number 
of fiare events to occur. They are scattered across the disk 
and we model this distribution by drawing the number of 
flares from Poisson's distribution around A^mean- To this 
aim we assume all flares last for 

Therefore, each simulated observation will show clear 
dispersion of the total luminosity around the assumed 
mean value, Lx- We can calculate this dispersion analyti- 
cally, deriving the mean number of flares during observa- 
tion A^mean from Eq. . Luminosity of a single flare is 
specified by Eq. (0) and combined with Eqs. I@J and Hll() . 
assuming random fluctuations around the mean value of 
flare number, as well as random variations in the distri- 
bution of flares and spots over the disk surface. 

The time-integrated luminosity is given by 



L{t) dt 



N 



(15) 



't=o 



where A is a random number drawn from the Poisson dis- 
tribution around A^meam and Ri and 4>i are random num- 
bers, which can be derived from the probability distribu- 
tion lO. 

The mean (expected) value of this quantity is defined 
in the same way as in Eq. (jHJ. This expression simplifies 
to 



L{t)dt ) = A, 



off 



-Rou 



mean 

.lif, , p(i?i)Li(i?i)di?i, (16) 

-Rin 



because t^^^ is constant, in accordance with the assump- 
tions. It is also evident that 




L{t)dt) -LxTobs. 



Variance of this quantity is defined analogously. 



(17) 




L{t)dt 



t=o 



it=o 



N=l 
N 



flout 



flin 



J2tMcLiiRu(bi) - LxTobs piRi)...piRN) 



i=l 




p{(j)i)...p{(j)N)dRi... dRiM d0i . . . d(t)N ■ 



(18) 



This integration can be carried out in a similar way as in 
Eq. JHJ, taking into account the relation 



J2 P{N;N^,,^)N^ = 7V„,ea„(iV„,ean + 1)- 



(19) 



The factor A^^ on the left side of this equation comes from 
the number of mixed terms, which arise from the expres- 
sion in square parenthesis. When the square of the first 
sum in square parenthesis is computed, the number of 
mixed terms is equal to N{N — l)/2. Therefore, the entire 
expression reduces to a single integral: 

/■flout 

= {t'igfN^,,,,J p{R,)LfdRi. (20) 

"'flin 

It is more convenient to use the normalized variance, i.e. 
the fractional variability amplitude i^var (see e.g. Vaughan 
et al. 2003) 

i^var = (21) 
^X ot)S 

which can be expressed as 



= 

var 



1 /;;7p(i?oi?di?i 



A^n 



jl^-'p{Ri)Li di?i 



(22) 



using Eqs. JH) and H20|l . In our case this expression is equal 

2 _ (/3rad- 2)2(1 -C2A--2)(1-C2) 1 



2C2(2/3,ad-2)(l-C'^"--2)2 



N„ 



(23) 
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The formula is more complex than the simple l/iVmcan 
dependence due to fact that different spots give different 
contributions to the total luminosity (parameter /^^ad de- 
scribes the flux energy distribution). It reduces to the fa- 
miliar case F^ar — V^mcan if ^11 flarcs have the same lumi- 
nosity (i.e. /3rad = 0). However, it shows strongly enhanced 
variability if flare luminosity scales with the flare radius 
proportionally to the dissipative disk flux (i.e. /Ji-ad = 3) 
since it this case the formula gives approximately 



Rn 



Rn 



1 



for i?out >> -Rin- 

We stress that Eq. (|23|l for the normalized variance 
does not include relativistic corrections. It can be applied 
to the case of a non-rotating Schwarzschild black hole with 
acceptable accuracy, but it turns out to be quite inaccu- 
rate for maximally rotating Kerr solution. 

Naturally, spots should only be generated between the 
adopted inner radius of the disk, i?in + Rx, and the outer 
radius, i?out — ^^x- In order to model a source with a 
prescribed luminosity, non-zero spot size has to be taken 
into account. A correction factor is obtained numerically 
through iterations. The spot radius Rx is constrained to 
be much smaller than i?outj and also it must be smaller 
than i?in. 

2.4. Generating an exemplary flare/spot distribution 

A specific distribution of spots can be generated from the 
adopted parameters of the model in four steps: 

(i) choice of the black-hole mass, mean source luminosity, 
duration of observation, and mean properties of the 
flare distribution, i.e. input model parameters M , Lx, 

-^obs 7 ^mcan ; -^in 7 Rout 5 -^0 7 Prad : and ^lifc , 

(ii) a starting guess of the spot radius Rx from Eq. 
and subsequent iterations of the radius, which allows 
to reproduce the required luminosity from the flare 
distribution between R^ + Rx and i?out — ^x; 

(iii) generation of the specific number of flares and corre- 
sponding spots, N, assuming the Poisson distribution 
with the mean TVmoan; 

(iv) generation of the position of each spot from the uni- 
form random distribution over the disk surface, again 
within the range (i?in -I- Rx, Rout — Rx)- 



The last operation means that the azimuthal angle 
a flare is determined as 



(jji = rnd X 2tt, (25) 
and the radius Ri is 
= [(-Rill + Rx)^ 

+ rnd X [(i?out - Rxf - (Rin + Rx)X^^ (26) 

where "rnd" is a random number between and 1. 

An example of the generated spot distribution is shown 
in Fig. n Each spot is localized at the disk surface at 



the moment of its generation. Spots occurring closer to 
the black hole are brighter than those which happen to 
be born farther out, because we assumed the power law 
scaling with /3rad > in Eq. 

2.5. Time-integrated properties of spots 

Spots are in Keplerian motion around the center. During 
the whole observation each spot moves along the (j)- 
direction, completing the distance 



(24) Acj^{r)^nK{r)t'^l. 



(27) 



Therefore, the integrated emission comes effectively from 
elongated belts rather than circular spots. Width of these 
belts is determined by radius of the flare and the length 
of arc which it has circumscribed, rA(f)(r). 

The belts are relevant for observational properties in- 
ferred from the model, because actual observations have, 
indeed, some finite (and non- negligible) duration, Tobs- If 
the integration time is long enough, some spots may com- 
plete one whole revolution, or even more full orbits, and 
so the length of the belt can be much longer than cir- 
cumference of a circle at the corresponding radius. This is 
particularly relevant for spots in the innermost part of the 
disk where motion is fast. Right panel of figure^shows an 
example of the resulting coverage of the disk surface with 
radiation belts. 

Computation of the time integrated spectrum is there- 
fore performed as the integration of local intrinsic emis- 
sivity of the belts. The locally emitted flux must be renor- 
malized in order to take into account that the emission 
now originates from the belt instead of a spot. Also, belt 
size is a function of radius Ri. We find 



-Fbclt — 



ttRI 



2i?xi?iA0(i?i) ^^^^ 

The spectral shape of the local spectrum is assumed to 
be the same for all spots (see Section 12.111 , and the re- 
flected spectrum is supplemented with the primary emis- 
sion in the proportion, which is determined by the radia- 
tion transfer computations fSect. m^i . The local emissiv- 
ity is assumed to be isotropic for both the primary flare 
emission and for the spot/reflected emission. Finally, ob- 
served time- integrated spectrum is subject to the relativis- 
tic corrections. 



of 2.6. Computations of general relativity effects 



In order to include general relativity effects in the Kerr 
spacetime, we have developed a new computational rou- 
tine which combines advantages of different approaches 
used in the past (Laor 1991; Karas, Vokrouhlicky & 
Polnarev 1992; Martocchia, Matt & Karas 2000). A re- 
design was desirable in order achieve sufficient resolution 
in both energy and time for the signal arising in multi- 
ple spots that are spread over the whole range of radii. 
The routine, ky, is flexible enough to allow easy modifi- 
cations of local emissivity profiles and it can be used as a 
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Fig. 1. Illustration of the spot distribution created by flares which illuminate the underlying disk. The plane of the disk 
is shown, as viewed from top, along the symmetry axis; (x, y) coordinates are scaled with the black hole gravitational 
radius, Rg. Left: a snapshot of an instantaneous spot distribution. Right: effective covering of the disk surface by 
radiating belts in a time integrated observation. Different realizations of the spot distribution have been explored in 
the paper (see Table EJ. 



rapid stand-alone code as well as linked to the standard 
XSPEC package (Arnaud 1996) for X-ray spectral analy- 
sis. In this paper we linked ky with the above-described 
computations of intrinsic emissivity of the spotted disk. 

The routine employs pre-stored data tables to look up 
and interpolate four relevant quantities. These are (i) the 
energy shift of photons, (ii) relative delay of their arrival 
time, (iii) magnification of radiation (lensing effect) in the 
Kerr metric, and (iv) the local emission angle (to account 
for directional anisotropy of the emission).^ Given the lo- 
cal (intrinsic) emissivity in the disk plane, the predicted 
(observed) spectrum can be integrated. In order to achieve 
higher accuracy than it was possible with previous rou- 
tines, the resolution in energy and time can be controlled 
and the grid covering the disk plane can be adjusted. Light 
rays are integrated in Kerr ingoing coordinates, so that 
the frame-dragging effect is tracked even very close to a 
fast rotating black hole (Misner, Thorne & Wheeler 1973). 
The equation of geodesic deviation is integrated separately 
to ensures precise determination of the flux magnification 
due to lensing near caustics. 



^ In fact, KY transfers all four Stokes parameters, which are 
necessary to compute light signal from a polarized source, 
which would be received in a detector equipped with a polar- 
izer. However, polarimetric information is superfluous for the 
present paper, and so we ignore it here. 



In general, the disk emissivity is characterized by radi- 
ation energy _Fem and the related photon fiux produced in 
the equatorial plane as function of emission energy Ecm, 



Fcni{E, 



cm; ' cnii 



(29) 



where -F'(rem) is total radiation flux emitted at the disk 
surface, ipi{Ecm) is the emissivity profile in frequency, 
ip2 (p) is the limb-darkening law (/i denotes the cosine an- 
gle between a ray and direction normal to the disk in the 
local co-rotating frame). In order to evaluate the terms 
in Eq. (|29() . one can employ analytic formulae or tabular 
data (both forms are used in this paper). In the case of a 
variable source in Eq. time-delay between light rays 
is taken into account at the point of their interception in a 
detector. Ray-tracing was performed in the Kerr geometry 
with redshift function g being given by 



and 



-g'^'^Cm + g^'^'TlrPr/Pt : 

Pan°' 



PaV 



(30) 



(31) 



Here, p" and rj" denote, respectively, the four-momenta 
of the photon and of the emitting material in the disk, 
and analogously p" and f)" that of the photon and of the 
observer at r — > oo; n" denotes a unit space-like vector 
perpendicular to the disk surface. For further details, see 
Dovciak et al. (2003). 
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As mentioned above, the influence of relativistic cor- 
rections on the local spectrum can be parameterized by 
the angular momentum of the black hole, a, and the in- 
clination angle of an observer, i. On the other hand, by 
scaling lengths with i?g and time intervals with i?g/c, one 
conveniently ensures that graphs of predicted spectra do 
not explicitly show dependence on the black hole mass. 

Relativistic effects, particularly in the case approach- 
ing a maximally rotating Kerr solution, significantly affect 
the shape of the radiation spectrum (especially at large 
inclination angle) . Also the magnitude of observed lumi- 
nosity as well as the observed variance turn out to be more 
sensitive to inclination if relativistic effects are taken into 
account. 

2. 7. Fractional variability amplitude 

We computed several radiation spectra representing dif- 
ferent realizations of the same mean spot distribution. In 
this way we simulated independent time-integrated obser- 
vations of a given source. We generate each distribution 
separately of the other, which means that we assume the 
lapse of time among simulated observations to be larger 
than the life-time of flares. 

Fractional variability amplitude was then calculated in 
standard way (see e.g. Edelson et al. 2002), as the energy- 
dependent variance of radiation flux normalized by the 
mean flux at a given energy. The variance was determined 
numerically, including relativistic corrections, from a finite 
number of realizations, K. Then the relative error of the 
normalized dispersion a was estimated as 

This dispersion error is purely due to intrinsic statistics 
properties and does not contain any additional finite count 
rate term usually present in the observational data. This 
specific formula is appropriate for our simulation because 
we draw the number of flares from Poisson's distribution. 
Models aimed at reproducing the entire power spectrum 
of AGN require more sophisticated approach. 

3. Results 

3.1. Locally emitted spectrum 

Radiative transfer computations were performed for flare 
parameters similar to those used by Ballantyne, Ross & 
Fabian (2001). We adopt the black hole mass of lO^M©. 
The irradiating X-ray flux, Fy^ — 10^^ erg cm^^ s~^, 
extends from 1 eV to 100 keV in the form of a power 
law (photon index F — 1.9). Iterations between the ver- 
tical disk structure (in hydrostatic equilibrium) and ra- 
diative transfer are performed till convergence is achieved 
(Rozanska et al. 2002). We choose the dimensionless ac- 
cretion rate of rri = lO^'^, (m is scaled to the Eddington 
rate, 2.2 x lO^^M/M© g s"^), and radius r = 18i?g as 
typical values where local computations are carried out. 



Hence, the ratio of the X-ray incident flux to the disk 
flux due to the internal dissipation for assumed incident 
flux and accretion rate is equal to 144. Disk thickness is de- 
termined self-consistently by proceeding from the disk sur- 
face down to the equatorial plane and adopting the diffu- 
sion approximation deeply below the disk atmosphere (for 
Thomson optical depth ^ 4) . The intrinsic spectrum of a 
spot (i.e. the reflection component) is shown in Figure 13 
Notice that the disk surface is highly ionized, and so the 
reflection is efficient also at low energies. 

Emission line properties from local spectrum 

Number of soft X-ray emission lines are visible in spec- 
trum, despite that the adopted irradiation is very strong. 
Since the soft X-ray emission lines have started to be 
seen in recent high-quality XMM and Chandra data for 
AGN (Kaastra et al. 2002, Rozanska et al. 2004), we list 
in Table Q the equivalent widths of the strongest lines, 
as measured with respect to the reflected continuum. The 
equivalent widths are determined with respect both to the 
reflected continuum and to the total continuum, i.e. pri- 
mary plus reflected. 

Sulphur and silicate lines are the strongest. Several 
other lines from highly ionized species, such as oxygen 
and carbon, are also well visible. 

An expanded region of the iron line (in linear scale) 
is shown in the right panel of Fig. El The iron line is 
strong and multi-component. This emission is mostly due 
to helium- like and hydrogen-like ions. The two-component 
structure is clearly visible, although ~ 6.7 keV component 
slightly dominates. Narrow-line components are accompa- 
nied by a broad shoulder due to Comptonization in the 
disk surface layers. The iron edge is quite deep, as usual in 
the case of highly ionized medium, and complex, as we see 
actually a two separate edges due to Fexxv and Fexxvi, 
at 8.83 and 9.28 keV correspondingly. These features are 
also considerably smeared by Compton scattering. 

3.2. General relativity effects in predicted spectra 

Relativistically smeared reflection from a disk uniformly 
covered with spots 

Thanks to its featureless power-law character, relativistic 
effects do not change the spectral shape of primary con- 
tinuum emission. However, spectral lines are still subject 
to the well-known relativistic smearing (e.g. Fabian et al. 
1989; Laor 1991). It is interesting to see if the smearing 
of a realistic spectrum is distinguishable from the case 
of separate monochromatic lines which may be superim- 
posed on broad-band continuum. To examine the effect 
of the smearing we consider in this section the disk sur- 
face, which is covered uniformly by a spot-like emission. 
The local emissivity is given by the reflection component, 
as specified in Sect. 13.11 The emissivity is assumed to de- 
crease with radius (/3rad = 3). 
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Fig. 2. Local reflection spectrum of an individual hot spot, which arises due to irradiation of the disk surface by a flare. 
Left: broad energy range is shown (log-log scale) with the intrinsic spectral resolution E/5E ~ 30. Dotted line shows 
the incident spectrum. Right: an expanded part of the same spectrum is plotted around the iron line energy (linear 
scale). Here, a higher intrinsic spectral resolution was adopted, E/SE = 100. In both panels, common parameters are 
M = IO^Mq, rh = 0.001, R = 9i?schw, = 10^^ erg s'^ cm'^, F = 1.9 (see Sect.Ofor details). 



There is a clear distinction between the relativistic 
smearing of a monochromatic line and the shape of rel- 
ativistically smeared realistic reflection component. From 
Figure 121 we see that the sharp blue edge, which is char- 
acteristic for a monochromatic line, is replaced by much 
smoother feature. This is due to the presence of double- 
peaked and Compton-broadened iron line in the realistic 
spectrum, as well as due to the contribution of the com- 
plex and Compton-smeared iron edge. 



Weak emission lines of other elements are also notice- 
able in smeared spectra. Therefore it is not surprising that 
traces of these relativistically broadened emission lines are 
noticeable in the observational data (e.g. Rozahska et al. 
2003 for Ton S180; Kaastra et al. 2002 for NGC 5548). 
Soft X-ray lines as strong as those claimed to be detected 
in the data by Branduardi-Raymond (2001) and Mason et 
al. (2003) are not expected within the frame of our model. 



This result clearly supports the argument (Zycki et 
al. 1997; see also Bao et al. 1998; Young et al. 1998, 
Martocchia et al. 2000; Gondoin et al. 2002; Ballantyne et 
al. 2003) that proper interpretation of the observational 
data requires consideration of the relativistic broadening 
of the entire reflected component, i.e. the reflected contin- 
uum together with the iron line, instead of treating the 
two components separately. 



Dependence of the luminosity on the inclination angle 

The importance of the relativistic effects for variability 
models was discussed in several papers (e.g. Abramowicz 
& Bao 1994; Xiong, Wiita & Bao 2000). Its complexity 
is mostly due to the effect of light bending and Doppler 
boosting. Relativistic effects play a crucial role in the 
observed time-dependent behaviour, particularly for the 
maximally rotating Kerr case. Therefore, in this section 
we start by considering all effects in this extreme case. 

The radiation originating in the innermost part of the 
disk shows significantly different dependence on the in- 
clination angle than the emission of the outer parts. We 
illustrate this by assuming, as a special example, that the 
disk surface is radiating uniformly in the local frame co- 
rotating with the disk. No separate spots and no radial 
dependence of luminosity are present. We consider two 
examples of such spectra: (A) the local emission is a pure 
power law with an energy index a = 0.9, (B) the local 
emission is a pure reflection component as described in 
Sect. 13.11 In both cases, we consider the emission coming 
from three detached regions: (1) an inner ring, r < 7i?g, 
(2) an intermediate ring, 7i?g < r < 20i?g, and (3) an 
outer ring, 20i?g < r < 50i?g. We assume no limb darken- 
ing (locally isotropic emission). 

In Figure0|we consider a strictly power-law emissivity 
profile and we show the corresponding behaviour of radia- 
tion flux at 6 keV on the inclination angle. We notice that 
the emission from the outer ring roughly follows purely 
geometrical effect which (trivial in the flat space), i.e. the 
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Table 1. Equivalent widths of most intensive emission 
lines from the local spot model (this paper), measured 

(i) with respect to the reflected continuum, EWrofl, and 

(ii) to the sum of the incident plus reflected continuum, 
EWtotai- Only lines with EWs higher than 1 eV are listed. 



Two forbidden Fexxv lines are transitions 2p ^P" — » Is^ 
and 2p — > Is^ ^So, Fexxv resonance line is the transition 
2p — > Is^ ^So- Two soft X-ray Fexxiv lines are transitions 



3d + 3s 
correspondingly. 



2p and 4d + 4s 2p ^P°, 



decrease proportionally to cos i, where i is the inclination 
angle of observation. However, the emission from the in- 
nermost part is more complicated: moving the observer 
from the disk axis towards its plane (i.e. i increasing), the 
observed flux first brightens, peaks at some intermediate 
inclination, and finally it decreases when the disk is in- 
clined very strongly. This behaviour is due to competition 
between the relativistic boosting and reduction of the pro- 
jected radiating area. The plot does not depend on energy 
at which the monochromatic luminosity is measured, be- 
cause the shape of the observed spectrum does not depend 
on inclination in the absence of spectral features. 

In the case of reflection component the results are prac- 
tically identical outside the region of intense iron line. On 
the other hand, there is more complex behaviour within 
the line range, as can be seen from Figure 13 

The assumption of radius-independent emissivity (be- 
tween i?in and i?out) and zero emission beyond the outer 
edge is rather artificial. Therefore we also show the re- 
sults for the disk whose radial emissivity normalization 
decays with the radius as (xr~^ (r < 50i?g). In Figure 
observed luminosity is plotted as function of inclination. 
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Fig. 3. The reflection spectrum from the spotted disk. The 
local intrinsic emissivity decreases (/3rad = 3) as radius 
increases up to i?out = 50i?g. The case of a Schwarzschild 
black hole, a = 0, and moderate inclination, i = 30° (panel 
a), can be compared with the case a = 0,i = 60° (panel b), 
rapidly rotating case a = 0.998, i = 30° (panel c), and a = 
0.998, i = 60°. Dotted lines show the case of relativistically 
smeared monochromatic line at 6.4 keV superimposed on 
a flat continuum. 



We see that the dependence is now even stronger than for 
the emission from the inner ring in Fig. 0] because with 
the adopted steep emissivity law most of the radiation is 
generated below 3i?g. Decline in luminosity is observed 
only at very large inclination, i.e. close to the disk equato- 
rial plane {i > 85°, probably unrealistic for real objects). 
Additionally, at such inclinations very close to 90° multi- 
ple images appear (Zakharov & Repin 2003). 
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Fig. 4. The dependence of the luminosity of the uniformly 
radiating disk belts on the inclination angle for the maxi- 
mally rotating Kerr black hole a = 0.998. Power law spec- 
trum and no radial dependence of the intensity were as- 
sumed. 

This dependency on inclination angle is only seemingly 
in contradiction with several other examples of such trend 
published in the literature (see e.g. Gierlihski et al. 2001, 
Ebisawa et al. 2001, Bhattacharyya et al. 2001, Ebisawa 
et al. 2001). There are three reasons why our results are 
different. 

First, the emissivity law is usually assumed to be less 
concentrated towards the center. When the emission from 
the multi-black-body Keplerian disk around a low mass 
object is calculated (like in all papers mentioned above) 
the emissivity goes down to zero at the marginally stable 
orbit due to the assumed boundary conditions. This is in 
accordance with the standard thin disk model. However, 
dissipation in the hot corona does not need to follow this 
law. Also the presence of the large scale magnetic field may 
eventually lead to more dissipation close to the marginally 
stable orbit (or even below it) than predicted in models 
based on zero-torque condition (sec Afshordi & Paczyhski 
2003 for the discussion). Models of the broad iron line 
actually indicate large concentration of the generated ra- 
diation flux towards the center (e.g. Wilms et al. 2001, 
Martocchia et al. 2002). 

Second, many computations have assumed some kind 
of limb darkening. For illustration we include in Figure |31 
the same limb-darkening profile as adopted by Laor 
(1991), / oc 1 -I- 2.06//. We see that in this case, under 
the assumed radial emissivity law, the fiux still rises with 
the inclination but more slowly. However, limb darken- 
ing characteristic for non-illuminated stellar atmospheres 
is actually replaced with limb brightening in the case of 



strongly irradiated disk atmosphere. The angle-dependent 
computations of the Compton reflection component show 
this effect clearly. The increase of emissivity with the in- 
clination strongly depends on the ionization state of the 
reflecting material. For the incident angle of 60° of the in- 
coming beam the ratio of local intensity at 90° to that at 
0° at '--^ 6 keV is equal to 9 (for weakly ionized medium) 
and 5 (for strongly ionized medium; cf. Zycki & Czerny 
1994; George & Fabian 1991; Martocchia et al. 2000). 
Therefore, in FigureElwe show also an example of the fltix 
dependence with limb brightening /(/i) oc 5 — 4/i roughly 
representing realistic reflection by the strongly ionized 
medium. In this case the flux rises almost by 50% with 
the change of inclination angle from 0° to 30°. However, 
in the total spectrum the primary emission dominates over 
the reflected component. The angular distribution of this 
emission is not known so in further considerations we as- 
sume isotropic emission. 

The dependence of the observed radiation on the in- 
clination angle has two important consequences. First, we 
recognize the tendency of the spectral shape on the incli- 
nation angle, as already illustrated in FigureOl Second, at 
low inclination angles the emission from innermost part 
of the disk is under-represented in comparison with the 
emission from outer parts, while at higher inclinations the 
opposite tendency appears. Since in the model outlined 
in our paper most of the variability is due to rare but 
luminous flares from the inner region, we flnd strong de- 
pendence of the variability properties on the inclination 
angle. We will discuss this in more detail in Sect. 13.31 

3.3. Fractional variability amplitude 

Parameters of the flare distribution 

The model depends on a number of parameters but most 
of them are in fact well constrained by observations, and 
hence the actual freedom is not so large. 

The mean number of flares at any moment, rimcan, can- 
not be less than ~ 10 because observed lightcurves do not 
show strong dips even on short timescales. Neither it can 
be very large, because in that case the variability ampli- 
tude would be far too small. In order to quantify variabil- 
ity, we use standard rms as its measure. Typical normal- 
ized rms of AGN is known to be of order of 0.3 (see e.g. 
Abrassart & Czerny 2000; Uttley et al. 2002; Markowitz 
et al. 2003). The parameter /3rad is expected to be ~ 3 
on the basis of the theoretical arguments as well as suc- 
cessful fits of some iron Ka profiles in such sources, in- 
cluding MCG-6-15-30 (e.g. Iwasawa et al. 1996, Sulentic 
et al. 1998, Nandra et al. 1997, Guainazzi et al. 1999), 
although more complex distributions were also advocated 
(see Merloni & Fabian 2003 and the references therein). 

In Table [3 we give several combinations of parameters 
that we considered together with the properties of the cho- 
sen flare distribution. In most cases, we fix the mass of a 
black hole at a representative value lO^M©, for which the 
local computations of the disc structure and radiation re- 
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Table 2. Exemplary flare distribution properties for the black hole/disk systeni seen at a moderate inclination, i = 30°. 
Parameters above the double line are assumed while the parameters below the double line are calculated as described 
in Sect.lO 



Parameter Model 







A 


B 


C 


D 


E 


F 


G 


I 


J 


K 


L 


P 


Q 




M 


10* 


10* 


10* 


10* 


10* 


10* 


10* 


10* 


10* 


10* 


10^ 


10^ 


10^ 




Rin 


6 


6 


6 


1.2 


6 


6 


6 


6 


6 


6 


1.2 


1.2 


1.2 




Rout 


50 


50 


50 


50 


100 


50 


50 


50 


50 


50 


50 


50 


50 


B 


/3rad 


3 


3 


3 


3 


3 


2.5 


3 


3 


4 


3 


3 


3 


3 


3 


Lx 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.1 


0.5 


0.5 


0.5 


< 


^mcan 


30 


30 


10 


100 


30 


100 


30 


100 


100 


30 


3000 


1000 


300 




ilifc 


10^5 


10^ 


10^ 


10'^ 


10^5 


10^5 


10^ 


10'^ 


10^5 


10^ 


10^5 


10^ 


10^ 




Tobs 


lO'' 


10^ 


10^ 


10^ 


lO'' 


lO'' 


10* 


10^ 


10^ 


10^ 


10^ 


10^ 


10^ 




Rl 


2.65 


2.65 


4.59 


0.62 


5.16 


1.57 


2.65 


1.45 


1.16 


0.84 


0.24 


0.43 


0.80 






3.17 


3.17 


6.11 


0.79 


7.42 


1.67 


3.17 


1.68 


1.32 


0.89 


0.28 


0.52 


1.10 


a 


mean 


60 


3030 


20 


200 


60 


200 


33 


200 


200 


60 


3030 


1010 


300 


B 


Cmcan 


0.12 


0.12 


0.15 


0.025 


0.16 


0.11 


0.12 


0.10 


0.07 


0.010 


0.09 


0.11 


0.14 


+J 


pa 

var 


43 


6.0 


74 


107 


81 


18 


58 


24 


34 


43 


27 


47 


87 


Q 


var 


28 


4.0 


38 


64 


36 


15 


38 


19 


26 


37 


22 


33 


45 




var 


26 


3.8 


39 


27 


36 


13 


36 


17 


25 


32 


7.0 


15 


28 



M is the black hole mass in PIq; Rin, Rout and R^ are the disk and the spot characteristic radii in Rg = GM/(? (radius _Rx of 
the spot orbit follows from Eq. mn . while i?x includes correction on finite size of spots); luminosity Lx is in units of 10** erg 
s~^ cm~^; time intervals tiife, Toha are in seconds; variances _Fvar are in percent; (Kfar follows from Eq. I)23|l . Kjlr includes the 
effect of the spot finite size and JVar was obtained numerically at 1 keV with relativistic corrections and finite size of the spots 
taken into account). 



processing were performed. The inclination angle of 30° 
was taken as a representative value (the view of the nu- 
cleus is frequently obscured at high inclinations by the 
presence of a molecular or dusty torus, so the average in- 
clination in type 1 objects must be significantly lower than 
60°). The opening angle of the dusty torus is statistically 
constrained by the ratio of type 2 to type 1 objects. It 
was estimated to be most probably between 2 and 4 (see 
Krolik 1999; Veron-Cetty & Veron 2000). Therefore, the 
torus opening angle is of order of ~ 40°, a typical type 
1 object is seen at ~ 30°, and a typical type 2 object at 
- 70°. 

Dependence of i<Var on energy 

An example of the X-ray spectra from two random real- 
izations of model A is shown in Fig.El The intrinsic spec- 
trum of spots/flares was assumed to be a combination of 
the intrinsic flare emission and reflected/reprocessed spot 
emission. The observed spectral features are now weaker 
than the example in Fig. [5] where only the spot emission 
was taken into account. The two spectra of statistically 
identical distributions of spots differ by normalization, but 
otherwise their form is very similar. Variations are slightly 
enhanced at energy where the strong iron line occurs. Tails 
of the line vary more significantly because this emission 
comes from several most luminous spots at the innermost 
part of the disk, and so random variations in the position 
of those few spots are relatively frequent and important. 



This effect is best seen in plots where the energy- 
dependent fractional variability amplitude is calculated 
from many random realizations. Such a plot (for model 
A) is shown in Fig. [7| 

In the case of a Schwarzschild black hole, variations 
are almost independent of energy. There is only marginal 
enhancement of variability in the red wing of the iron line, 
around 5 keV. There is also a small dip in the plot right at 
~ 7 keV. However, the overall energy dependence is not 
strong, and so the value of the dispersion at 1 keV is quite 
representative for the entire spectrum (Table 

The energy dependence of Fvar is stronger in the case 
of emission from the disk surrounding a fast-rotating Kerr 
black hole. We show the corresponding plot for model D 
in the lower panel of Fig. |3 Here we adopted a higher 
value of the mean number of clouds, because otherwise 
the obtained variance would be unrealistically high in the 
context of AGN. We see in this case strong variations in 
the region of the red wing of the iron line, around 5 keV, 
and again a small dip in the plot right at ~ 7 keV. Still, the 
overall variability is well characterized by the dispersion 
at 1 keV, as given in the Tab. El 

The fractional variability amplitude must be in this 
case calculated with the relativistic corrections. The ef- 
fect depends also on the inclination angle. An example for 
the model D but computed this time at much higher in- 
clination angle of 70° is also shown in Fig. [7| At such a 
high inclination, line wings are still broader than at 30°, 
influencing the entire energy range between 4 and 8 keV. 
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Fig. 5. The dependence of the luminosity of uniformly ra- 
diating disk on the inclination angle with isotropic local 
emission, limb brightening (/(/i) oc 5 — 4^), and limb dark- 
ening (/(/i) oc l-|-2.06/i). The radial emissivity law oc r^^ 
was assumed together with a power-law spectral shape. 



Dependence of Fvar on the inclination angle 

We have noticed already in Fig. the strong dependence 
of the overall variance on inclination angle, which occurs 
in the case of rapidly rotating Kerr solution. In order to 
examine this trend in more detail, we plot in Fig. [S] the 
variance measured at 1 keV as a function of the inclination 
angle. 

There is a noticeable difference between the variability 
level of the source observed at i = 0° versus the case of 
i — 30° (the value, which is typically expected for broad- 
line galaxies on the basis of AGN unification scheme). 
The fractional variability amplitude is enhanced by 47%. 
Further increase, by a factor of 3.1-4.6 in the variability 
level, is expected for objects which are observed at i ~ 60- 
77°. Such inclinations are typical for Seyfert 2 galaxies. In 
those objects, in the soft X-ray band, we do not see the 
direct emission from the nucleus. Instead, in some of them 
the direct emission is visible in hard X-rays (Done et al. 
2003), which opens a way to test our prediction regarding 
the expected level of variability. 

Several author recognized importance of general rela- 
tivity effects for the observed dependence of the variability 
properties on the inclination. Abramowicz & Bao (1994) 
noticed the flattening of the power spectrum slope with 
an increase of the inclination angle. Their approximate re- 
sults were in general confirmed by Xiong et al. (2000) for 
Schwarzschild geometry. Karas (1997) discussed Fourier- 
phase analysis of the spotted disk model. 



Fig. 6. Two random realizations of the same spot distri- 
bution (thin lines) and the original local emission (thick 
line histogram) for model A; the ratio of the two random 
realizations are shown below. The local input spectrum is 
the sum of the primary (flare) emission and the reflected 
(spot) component, without any limb-darkening. 




E [keV] 



Fig. 7. Fractional variability amplitude was calculated as 
a function of energy. The model A (i = 30°; upper panel), 
the model D (i = 30°; middle panel), and the model D 
at high inclination observed (i = 70°; bottom panel). 
Duration of a single observation 10^ s, number of expo- 
sures 300. 
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Fig. 8. The dependence of the normaUzed dispersion mea- 
sured at 1 keV on the inclination angle of observation 
for model A (upper panel, continuous line), model J (up- 
per panel, dashed line) and model L (lower panel, middle 
curve). Two other curves in the lower panel show the re- 
sult of increasing the mean number of flares to 3000 and 
decreasing the mean number of flares to 300. Number of 
exposures was 100. 



The quantitative effect of the inclination depends on 
details of the assumed variability model. We show this 
dependence for model J, which is characterized by steeper 
emissivity {Prad = 4) than what was adopted in other 
cases. We see that the rise of the fractional variability 
amplitude with the inclination is faster than for /3rad = 3 
because more radiation comes from the innermost region 
where the relativistic effects are important. 



The importance of the limb brightening 

As mentioned in Sect. 13.21 limb brightening is expected 
for the spot emission. We found, however, that it has 
no strong effect on the overall dependence of the frac- 
tional variability on the inclination angle (the increase of 
the variance with the inclination angle is still present, al- 
though it is slightly reduced) . Limb-brightening of the pri- 
mary emission (in the form adopted in Fig. results in 
further flattening of the plots. For the model P, the ratio 
of the variance at 60° to the variance at 1° is equal to 
3.1 while for the same model with limb brightening of the 
spot emission is equal to 2.6, and for the limb brightening 
of both the flare and spot emission is equal to 2.2. 



The value of the fractional variability amplitude given by 
Eq. H23|l is independent both from the energy and incli- 
nation angle. In TableElwe give these analytical values as 
f^a^j.. We also calculate semi- analytical values F^^,. taking 
into account the finite size of the spot but neglecting rel- 
ativistic corrections; a new radius is calculated iteratively 
and inserted into Eq. with ( = (i?in + Rx)/{Rout — 
i?x). 

We see that the analytical formula for F^^,. gives a 
rather good approximation of the exact value in the model 
K, which involves the Schwarzschild solution and rel- 
atively small size of spots (less than li?g). For other 
Schwarzschild models with larger size of spot, the sim- 
plest approximation is not satisfactory any more, but i^^^j. 
still provides a reasonable approximation. Relativistic ef- 
fects are important but the formulae neglecting them can 
be used for qualitative estimates. 

For Kerr models, the analytical value given by Eq. H23|l 
over-predicts the variability level at low inclination angle 
by a factor up to 3, but it represents roughly the correct 
values at large inclination angles. In this case, analyti- 
cal estimates give large errors, and so accurate numerical 
energy-dependent and inclination-dependent results with 
general relativity effects are needed. 

Trends with assumed model parameters 

We see that short duration of the flare (below a few 
xlO^ s) leads to too small normalized dispersion in X- 
ray lightcurves for a typical long observation lasting about 
1 day. It does not mean that short-lasting flares are ruled 
out, but such flares clearly cannot be responsible for vari- 
ability in the timescales of days or more. 

Increasing the assumed radius of the disk increases also 
the variance because, in our model, larger -Rout translates 
into a smaller number of flares generated in the innermost 
part of the disk where flares are the most energetic. This is 
clearly connected with our assumption of the uniform cov- 
ering of the disk and the power-law radial decay of flare 
luminosity, as given by /3rad- Equivalently, we can mod- 
ify the dispersion by changing the slope (irad- lower value 
(model F) reduces variations, while higher value increases 
it (model J). 



3.4. Special case of MCG-6-30-15 

This nearby (z = 0.007749; cf. NASA/IPAC Extragalactic 
Database) NLSl galaxy is one of the best studied ob- 
jects in the X-ray band. It is well-known for its extremely 
broad, relativistically smeared iron line (Tanaka et al. 
1995, Fabian et al. 1995, Iwasawa et al. 1996, Fabian & 
Vaughan 2003 and the references therein). The slope of 
the power law component in this source (photon index 
r = 1.87 ± 0.01; Wilms et al. 2001) is not significantly 
different from the value F = 1.9 adopted in our radiation 
transfer computations. 
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Energy dependent fractional variability amplitude was 
studied for this source by Matsumoto et al. (2003) on 
the basis of ASCA data. They showed that the overall 2- 
10 keV amplitude of ^ 18-20% is enhanced in the region 
of red wing of the iron line (modeled as a broad feature 
around 5.4 keV) up to 40% (but with large error), and 
somewhat less at its blue wing (around 6.4 keV). Detailed 
analysis was performed by Markowitz et al. (2003) on 
the basis of several RXTE observations. They found that 
the overall variability decreases with energy, going from 
~ 20% at 4 keV down to - 13-14% above 10 keV. There 
is also a slight relative enhancement at ~ 5 keV and at 
~ 8 keV, and a dip at ~ 6.5 keV. Similar trend was no- 
ticed by Fabian et al. (2002) in their long XMM observa- 
tions, but in this case the variability level was found to 
be slightly higher. We model the observed fractional vari- 
ability assuming parameters appropriate for this source. 

The mass of the black hole comes out rather small on 
the basis of X-ray power spectrum analysis: logM = 6.0 
follows from the normalization method (Czerny et al. 
2001), while logM = 5.7 has been derived from the high- 
frequency break model (Uttley et al. 2002). However, in 
view of uncertainties, the mass as large as logAf — 5.9 
is still compatible. No results have been obtained from 
reverberation studies for this source. The mass was esti- 
mated to be log A/ = 7.0 (Czerny et al. 2001) using the 
disk method. Reynolds (2000) also argues for this higher 
value on the basis of the bulge-black hole mass relation 
(Magorian et al. 1998) We thus adopt M = lO'^M© in our 
computations. 

X-ray luminosity varies significantly, but taking the 
typical 2-10 keV flux to be 4 x 10""erg s"^ cm"^ 
(Weaver et al. 2001) we can estimate the total X-ray flux 
(1 eV-100 keV) to be of order of 5 x 10'*'^ erg s~^ (assuming 
the Hubble constant Hq — 70 km Mpc~^). 

The bolomctric luminosity of this source can be esti- 
mated either from the 2-10 keV luminosity, taking into 
account the bolomctric correction of 27 (Padovani & 
Rafanelli 1998), or directly from the broad-band data cor- 
rected for extinction. The first method gives 1.6 x 10"*^ 
erg s""'^ while the second one gives 9.1 x 10^^ erg s~^ 
(Nikolajuk, in preparation). This means that for the 
adopted mass of the black hole the source radiates at 
~ 10% of the Eddington luminosity. 

We considered three values of the mean number of 
flares: 300, 1000, and 3000. The corresponding models are 
denoted L, P and Q in Table [21 We choose the duration 
of the observations to be 10'^ s, the typical duration of 
continuous observation in the data, but the result does 
not depend strongly on this value. The adopted life time 
of flares, 10^ s, is chosen to reproduce the knee in the 
power spectrum of this source, as determined by Uttley et 
al. (2002). Therefore, we model the day-to-day variability 
(and the result docs not apply to longer timescales). 

Comparison of the observed fvar in short (day-to-day) 
timescales (see Figure 5 of Markowitz et al. 2003) with 
theoretical level of variability indicates that the most ap- 
propriate mean number of flares is between 300 and 1000. 



There is evident similarity in the energy dependence of 
i^var between observed data and the theoretical diagram: 
variability in the region of the iron line red wing (~ 5 keV) 
is enhanced, and a dip occurs at about 7 keV. However, 
the dependence on energy is generally much weaker in the 
model than in the data. Also, the model does not show 
the observed strong rise of the i^var towards low energies. 

This presence of strong variation at ~ 3-4 keV in the 
data and the lack of such variations in the model may indi- 
cate that indeed (Inoue & Matsumoto 2003), this energy 
band is still affected by varying conditions in the warm 
absorber. 

4. Discussion 

The flare / spot model is an attractive explanation of the X- 
ray spectra and variability of AGN. In this scenario X-ray 
emission is generated both in hot magnetic loops above 
an accretion disks and in the bright spots created under 
the loops by strong irradiation. In the present paper we 
tested this model by analyzing the predicted fractional 
variability amplitude, i^var- 

We derive simple analytical formulae which allow to 
estimate the level of variability from the assumed mean 
number of flares, flare duration, integration time of a single 
observation, the ratio of inner to outer disk radius and 
the radial dependence of the flare luminosity. This energy- 
independent and inclination-independent formula roughly 
applies to the case of a non-rotating black hole and small 
size of the spots. The formula is based on assumption of 
the uniformly covered disk surface but it can be readily 
generalized to another kind of distribution. 

Larger spots and/or fast rotating black hole, with inner 
disk radius close to the marginally stable orbit require 
numerical approach, and we show both the dependence of 
the resulting _Fvar on energy and inclination angle. 

This leads to a firm prediction of our model which can 
be used to test the basic scenario. The model shows that 
if the disk extends close enough to the black hole, the 
general relativity effects lead to significant enhancement 
of the variability at large inclination angle of observation. 
Therefore, Seyfert 2 galaxies, if intrinsically identical to 
Seyfert 1 objects but viewed at a larger angle, should ex- 
hibit statistically higher Fvm- when measured at the same 
energy. A factor of 3 difference between Seyfert 2 and 
Seyfert 1 galaxies is expected for a Kerr black hole and 
a factor of 1.4 for a Schwarzschild black hole. 

The observational evidence of any trends in this di- 
rection is scarce at present, but it may support or reject 
our view more firmly in future. However, some values can 
be given already now. Mean variability level of Seyfert 1 
objects is ^ 18% on the short timescales of 1 day in the 
sample of Markowitz et al. (2003). Seyfert 2 galaxy NGC 
4945 is seen through the torus and it displays variability at 
the level of ^ 40% in the hard X-ray band. As another ex- 
ample - NGC 7582 (Mihara et al. 2000), has revealed the 
normalized variability amplitude in the hard X-ray band 
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to be about ^ 30%. Studies of more Seyfert 2 objects are 
clearly needed. 

The overall intra-day variability level of Seyfert galax- 
ies is well explained by the flare model if the mean number 
of flares is of the order of 30-100 for assumed non-rotating 
black hole, or 300-1000 for a fast rotating black hole. The 
number of flares requested is much larger than the usual 
expectation of ~ 10 flares. Such a small number of flares 
is predicted when an assumption is made that all flares 
have the same luminosity. In our more realistic model we 
allow for the flare luminosity to depend on the occurrence 
radius. When we assume that the flare luminosity scales 
with radius as cx r~^, most of the source X-ray luminosity 
comes from a few flares generated in the innermost part 
of the disk which enhances the variability. This trend is 
clearly seen from our analytical expression 1231 for the nor- 
malized variance, which, for Prad — 3, and i?out >> -Rin, 
reduces to Eq. Therefore, the total number of flares, 
-^moan Can Still bc quitc high for a source with a moderate 
variance. 

The energy dependence of -Fvar is generally weaker in 
the model than in data. Observed variations show trends 
with energy in 1-10 keV band. For example, in Markowitz 
et al. (2003) Akn 564 varies at the level 18-22%, IC 4329A 
at the level 11-14% and MCG-6-30-15 at 14-21%, depend- 
ing on the considered energy. In our models the trends in 
Fvar with energy never exceed 2%. 

This discrepancy can be possibly solved by relaxing 
several simplifications: 

— the local spectrum was computed in detail only at a 
single radius. In reality, the shape of the spot spectrum 
is expected to show significant trends with the disk ra- 
dius, as for example emphasized by Zycki & Rozahska 
(2001). A grid of spot spectra should be computed (pa- 
rameters space of the model is rather rich, and so the 
task is numerically very time consuming and we post- 
pone it for future work); 

— hydrostatic equilibrium was assumed to compute the 
irradiated disk structure. However, timescales of the 
flares and of restoring the hydrostatic equilibrium are 
roughly comparable (sec the discussion by Nayakshin 
& Kazanas 2002 and Colhn et al. 2003). Therefore, 
neither the assumption of hydrostatic equilibrium nor 
the assumption of unperturbed disk are satisfactory. 
Actual disk evolution should be followed; 

— we assumed a unique value of flare duration since 
we aimed at modeling the variations at the dominant 
timescale, e.g. at the knee of the power spectrum. A 
distribution of flare timescales, a coupling between the 
flare occurrence (avalanches), and an exact profile of 
an individual flare are needed if we want the model to 
reproduce the entire power spectrum (e.g. Lehto 1989, 
Abramowicz et al. 1991; Xiong et al. 2000, Merloni & 
Fabian 2001); 

— we neglected the possible effect of the variable warm 
absorber which may be important in the soft X- 



ray band for some sources, as argued by Inoue & 
Matsumoto (2003); 
— we neglected the possible contribution of the radiation 
reprocessed by some distant reflector, like an outer disk 
or dusty /molecular torus (e.g. Krolik et al. 1994). 

It is to be seen whether elimination of these assump- 
tions would lead to better agreement of the predicted en- 
ergy dependence of i^var with the data. Particularly diffi- 
cult seems to bc the explanation of both apparently lower 
observed variability in the iron line region as well as the 
enhancement of the variability towards low energies seen 
in Fig. 5 of Markowitz et al. (2003). 

We considered just a few special cases along this line. 

A change of the life time of a flare from radius- 
independent to scaled with Keplerian timescale, without 
a change of other parameters, resulted in a fractional 
variability amplitude even less dependent on the energy 
than previously. It is simply caused by the fact that in 
the case of such a scaling relatively more energy is dissi- 
pated in outer region, so the most relativistically broad- 
ened and variable inner region contributes less to the total 
lightcurve. The overall normalization of the -FVar depends 
on the proportionality constant between the life times and 
the Keplerian timescale (with other parameters flxed). 

If we assume that the overall radial dependence of the 
dissipation should not be modified, we can consider two 
representative examples. 

In the first case we assume that the life time of a fiare 
scales with the Keplerian timescale but the probability 
of a fiare to appear at a given radius scales inversely with 
the local Keplerian timescale. We have therefore less long- 
living flares in the outer region and more short-living flares 
in the inner region, with scaling of a single flare luminosity 
with radius unchanged. In this case again the fractional 
variability amplitude is less dependent on the energy than 
in our basic model. Having more flares localized in the 
inner region lead to reduction in the fluctuations in the 
relativistically smeared red wing. 

In the second case we again assume that the life time of 
a flare scales with the Keplerian timescale, we still adopt 
the uniform distribution of the flares across the disk sur- 
face but this time we assume that a flare luminosity de- 
creases with radius even more strongly (/3rad — 4.5 instead 
of usually adopted /3rad = 3) to compensate for an increase 
of the flare life time. Such a solution leads to slightly en- 
hanced dependence of the fractional variability amplitude 
on the energy but the effect is not strong. 

Complex dependence of the spectral shape of the re- 
flected component on the disk radius may introduce sig- 
nificant modification to the predicted energy dependence. 
Clear suggestion of what is needed can be seen from the 
recent analysis of MCG-6-30-15 by Vaughan & Fabian 
(2003). A constant component with complex energy de- 
pendence is apparently needed in order to formally model 
the fractional variability amplitude in this source (see their 
figs. 11 and 16, top panel). This component may perhaps 
be understood as a reflection component (above ~ 1 keV) 
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and a contribution from emission/scattering by some ex- 
tended medium (below ~ 1 keV). However, it is not clear 
whether there is any possibility to find a parameter range 
which would satisfy two 'opposing' trends seen in these 
data: we need strong reflection from distant disk region in 
order to reproduce the constant component but we need 
strong reflection from innermost region in order to explain 
the strongly relativistically broadened iron line profile also 
seen in these data. Attempts by Rozahska & Zycki (2001) 
and Ballantyne et al. (2003) were not successful. 

Flares are not the only possibility to explain the 
X-ray emission of AGN. Other scenarios include lamp- 
post (standing shock) model (e.g. Henri & Pelletier 1991, 
Malzac et al. 1998), model of gradual or rapid disk evap- 
oration and its replacement by the hot flow (e.g. Narayan 
& Yi 1994, Liu et al. 2002), possibly with an outflow (e.g. 
Blandford & Begelman 1999), and the cloud model (e.g. 
Collin et al. 1996, Karas et al. 2000). Further work, tak- 
ing into account variability issues, is needed to determine 
whether the flare/spot model is the most satisfactory. 
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